Mon. Not. R. Astron. Soc. 000,[T]-?? (2001) Printed 5 February 2008 (MN KTpX style file v2.2) 



Have we detected the most luminous ULX so far? 
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ABSTRACT 

We report the XMM-Newton detection of a moderately bright X-ray source (F .5_7 ~ 8.2 x 
10~ 14 erg cm -2 s -1 ) superimposed on the outer arms of the inactive spiral galaxy MCG- 
03-34-63 (z=0.0213). It is clearly offset from the nucleus (by about 19") but well within the 
D 2 5 ellipse of the galaxy, just along its bar axis. The field has also been observed with the 
Hubble Space Telescope (HST) enabling us to compute a lower limit of > 94 on the X-ray 
to optical flux ratio which, together with the X-ray spectrum of the source, argues against 
a background AGN. On the other hand, the detection of excess X-ray absorption and the 
lack of a bright optical counterpart argue against foreground contamination. Short-timescale 
variability is observed, ruling out the hypothesis of a particularly powerful supernova. If it is 
associated with the apparent host galaxy, the source is the most powerful Ultra-Luminous X- 
ray source (ULX) detected so far with a peak luminosity of ~ 1.35 x 10 41 erg s" 1 in the 0.5- 
7 keV band. If confirmed by future multi-wavelength observations, the inferred bolometric 
luminosity (~ 3 x 10 41 erg s~ 4 ) requires a rather extreme beaming factor (larger than 115) 
to accommodate accretion onto a stellar-mass black hole of 20 M® and the source could 
represent instead one of the best intermediate-mass black hole candidate so far. If beaming is 
excluded, the Eddington limit implies a mass of > 2300 Mq for the accreting compact object. 

Key words: X-rays: binaries - X-rays: individual: XMMU J132218. 3-164247 - galaxies: 
spiral - galaxies: individual: MCG-03-34-63 - black hole physics 



1 INTRODUCTION 

All known black holes belong to two families: stellar-mass black 
holes are seen in X-ray binaries, while super-massive ones are 
present in the centres of galaxy bulges sometimes revealing them- 
selves as Active Galactic Nuclei (AGN). While the former have 
masses up to ~ 2QMq (e.g. Fryer & Kalogera 2001), the lat- 
ter have masses in the range ~ 10 6 -10 9 Mq, the smaller-mass 
super-massive black hole to date being that in NGC 4395 with a 
mass of a few times 10 Mq (Peterson et al. 2005). Although it 
has long been thought that intermediate-mass black holes (IMBH) 
with masses ~ 10 2 -10 4 Mq may form in dense stellar clusters 
(e.g. Frank & Rees 1976; Portegies Zwart et al. 1999), there are no 
known IMBHs filling the mass-gap between the two known fam- 
ilies. If present, active IMBHs may reveal themselves as accreting 
X-ray sources exceeding by a large factor the Eddington luminos- 
ity of stellar-mass black holes (Lf o d< ^ = 2.6 x 10 39 erg s" 1 for 
a "maximal-mass" stellar-mass black hole of 20 Mq). 

Ultra-Luminous X-ray sources (ULX) are off-nuclear point- 
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like X-ray sources seen in other galaxies (than the Milky Way) 
with luminosities exceeding I^oMq (see e.g. Colbert & Mushotzky 
1999; Mushotzky 2004). Since the Eddington argument implies a 
lower limit of 20 Mq on the mass of the central object, ULXs 
are often regarded as IMBH-candidates (see e.g. Miller & Col- 
bert 2003; Fabbiano 2005). However, inferring a lower limit on 
the mass of an accreting compact object only from its bolomet- 
ric luminosity can lead to misleading results. This is because, if 
potential anisotropies of emission (e.g. beaming, see Reynolds et 
al. 1997; King et al. 2001) or accretion (e.g. radiation-driven in- 
homogeneous accretion, see Begelman 2002) are not taken into ac- 
count, the lower limit on the mass of the object may be severely 
over-estimated. Both the beaming and inhomogeneous accretion 
scenarios can be invoked to explain luminosities up to a few times 
10 40 erg s _1 with accretion on standard stellar-mass black holes. 
However, geometric beaming (e.g. a funnel geometry) and inhomo- 
geneous accretion can only provide an effective luminosity exceed- 
ing the Eddington limit by a factor ~23 (Madau 1988) and ~10, 
(Ruszkowski & Begelman 2003) respectively. Thus, for ULXs with 
luminosities exceeding 10 41 erg s" 1 , relativistic beaming seem the 
only option to avoid the presence of an IMBH. It is thus clear that 
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Figure 1. For all images, North is upwards and East to the Left. Top Left: a 3.5' X 3.5' DSS optical image in the region of IRAS 13197-1627 and MCG-03- 
34-63 is shown together with the detected X-ray sources (red circles). We also show the D25 ellipse of MCG-03-34-63 as a black ellipse. Top Right: the grey 
scale is the VLA 21 cm image in the region of of the galaxy MCG-03-34-63. Optical contours from the DSS are shown in blue. The red circle is centred on 
the XMM-Newton position of ULX 1 and the red cross is the optical centre of MCG-03-34-63. Bottom Left: a portion of the HST WFPC2 500s exposure 
(F606W filter) of MCG-03-34-63 is shown together with the XMM-Newton error-circle (2" in radius) centred on ULX 1, Bottom Right: expanded view of 
the HST image around the X-ray error-circle. The yellow box is centred on the brightest optical source within the XMM-Newton error-circle and represents 
the position of the most likely optical counterpart of ULX 1 . 



the best IMBH-candidates are ULXs in the high end of the lumi- 
nosity function with luminosity approaching (or possibly above) 
10 41 erg s _1 . The best IMBH-candidate so far is a ULX in the 
galaxy M 82. Its extreme luminosity (up to ~ 9 x 10 40 erg s -1 ), 
together with the detection of a 54 mHz QPO and an Fe Ka line 
(which both argue against a beamed source), make it difficult not 
to invoke an IMBH (Matsumoto et al. 2001; Kaaret et al. 2001; 
Strohmayer & Mushotzky 2003). 

Here we report the XMM-Newton detection of an X-ray 
source which is plausibly a ULX in the galaxy MCG-03-34-63. 
If it is at the distance of its apparent host galaxy, the bolometric 
luminosity of ~ 3 x 10 41 erg s _1 , exceeds by a factor > 115 the 
Eddington luminosity for a 20 Mq stellar-mass black hole, mak- 
ing it the most luminous ULX ever detected, and possibly one of 
the best IMBH-candidate so far. 



2 MCG-03-34-63 ULX-1 

XMM-Newton observed the active galaxy IRAS 13197-1627 



(z=0.0165) on 2005 January 24 for about 45 ks (Miniutti et al. 
2006). IRAS 13197-1627 is not isolated in the sky and the inac- 
tive spiral galaxy MCG-03-34-63 (z=0.0213) is only 1.8' away to- 
wards NW, in the field of view of all the XMM-Newton EPIC cam- 
eras. MCG-03-34-63 has been often confused with IRAS 13197- 
1627 (a.k.a. MCG-03-34-64) until 1996, where it was detected in 
the radio with a total flux of 93.4 mjy at 4.9 GHz (Colbert et al. 
1996). It is also a relatively luminous IR galaxy with I/4o-i22nm = 
2.1 x 10 10 1/0, most likely associated with star-formation (Surace, 
Sanders & Mazzarella 2004). 

In the top left panel of Fig. 1, we show the optical Digitized 
Sky Survey (DSS) image centred between the two galaxies. We 
also show as circles the X-ray sources detected in the field, and 
the D25 ellipse of MCG-03-34-63 (the dimension defined to be 
equal to 25 mag arcsec~ 2 ). Besides IRAS 13197-1627, the second 
strongest X-ray source is superimposed on the outer arms of the 
galaxy MCG-03-34-63 1 . Its X-ray position is displaced by about 



A second source (S2 in Fig. 1) is also detected. Its X-ray spectrum sug- 



© 2001 RAS, MNRAS 000,fT]-?? 



MCG-03-34-63 ULX 1 3 



19" in the SW direction with respect to the optical centre of the 
galaxy and its coordinates (J2000) are RA 13h22ml8.3s and DEC 
-16d42m47.7s. Its official name is thus XMMU J132218.3-164247 
and we shall call it ULX 1 here for brevity. The X-ray position of 
ULX 1 has been determined by running the source detection task 
edetect_chain and is corrected by cross-correlating 4 X-ray 
sources with the USNO-B 1 .0 optical catalogue (Monet et al. 2003). 
The la r.m.s. absolute astrometric accuracy of XMM-Newton 2 is 
about 2". Here, we make a conservative choice and we assume an 
error-circle of 2" in radius. The X-ray position will be soon known 
with much better accuracy thanks to an approved Chandra observa- 
tion. No optical point-like counterpart is seen in the DSS meaning 
that the source is fainter than ~21 in R and ~22.5 in B. 

We also analysed the Optical Monitor (OM) mosaic of the 
field images (UVW1 filter). The UV contours trace MCG-03-34- 
63, peak at about 17" NE of the galaxy centre, and no point- 
like source corresponds to ULX 1. Radio images (VLA 21 cm, 
1993 May) have been analysed as well and, besides the AGN 
IRAS 13197-1627, there is an excess within the D25 ellipse of 
MCG-03-34-63 which is however not consistent with ULX 1 nor 
exactly centred on the galaxy optical centre (see top right panel 
of Fig. 1). At 4.9 GHz, the nucleus of MCG-03-34-63 is clearly 
detected (Colbert et al. 1996), but there is no radio excess at the 
position of ULX 1 . 

The Hubble Space Telescope (HST) also observed the region 
in July 1994 for 500 s with the WFPC2 and the F606W filter ap- 
plied. A portion of the HST image is shown in the bottom left panel 
of Fig. 1 together with a white circle of 2" radius centred on the X- 
ray position of ULX 1 and representing the XMM-Newton error- 
circle. The HST astrometry has been corrected by using 5 common 
sources with the USNO-B 1.0 catalogue. It is interesting to notice 
that ULX 1 lies along the bar axis of MCG-03-34-63 where the 
interstellar medium may be stressed and where off-nuclear X-ray 
sources have been found in other galaxies (see e.g. in NGC 1672, 
Brandt et al. 1996). In the bottom right panel of Fig. 1, we show 
a zoom into the HST image. The brightest optical source within 
the XMM-Newton error-circle is marked with a yellow box and 
represents the most likely optical counterpart (if any) of ULX 1 . 
Another faint source is present in the error-circle and two are seen 
at its edges towards SE. 

In the X-rays, the field has been observed with ASCA and Bep- 
poSAX, but the moderate angular resolution of the two X-ray mis- 
sion detectors does not allow us to distinguish between the bright 
AGN IRAS 13197-1627 and ULX 1. IRAS 13197-1627 is de- 
tected by Swift/BAT (Markwardt et al. 2005). However, given the 
17' angular resolution nothing can be said on ULX 1. No pointed 
ROSAT observation exists and ULX 1 is too faint to be detected in 
the ROSAT All Sky Survey. 



3 X-RAY DATA ANALYSIS 

ULX 1 provides 457 background-corrected counts in the EPIC-pn 
camera, and about 300 in each of the MOS detectors in the 0.5- 
7 keV band. We extracted the X-ray light curve of ULX 1 from the 
MOS+pn detectors to search for short-timescale variability. Given 
the low count rate (about 0.025 Cts/s), we chose a relatively long 

gests it is a Compton-thin Seyfert 2 galaxy at z~0. 16 (Miniutti et al. in 
preparation). 

2 see Kirsh et al. 2006, XMM-Newton calibration document XMM-SOC- 
CAL-TN-0018 at|http://xmm.vilspa.esa.es| 
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Figure 2. We show the background subtracted light curve of ULX 1 ob- 
tained from the MOS+pn data in the 0.5-7 keV band (binsize of 2.9 ks). 
The Kolmogorov-Smirnov probability that the light curve is constant in 
5.7 X 10" 4 only. 
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Figure 3. The pn (black) and MOS (red and green) spectrum of ULX 1 
is shown together with the best-fitting absorbed power law model and the 
resulting residuals (see text for discussion). 

time bin of 2.9 ks in order to reduce the errors on the individual data 
points. The resulting background-subtracted light curve is shown 
in Fig. 2. A fit with a constant gives x 2 = 27.3 for 12 dof, show- 
ing the the source is variable at the 99.3 per cent confidence level. 
By applying the Kolmogorov-Smirnov test, the significance of the 
variability is increased to the 99.94 per cent level. 

We also extracted the X-ray spectra from both the pn and 
MOS cameras in the same 0.5-7 keV band. As a first attempt to 
describe the X-ray spectral shape of ULX 1, we consider an ab- 
sorbed power law model and we fix the absorbing column to the 
Galactic value, i.e. at 5.8 x 10 20 cm" 2 (Dickey & Lockman 1990). 
The joint fit to the pn and MOS data is unacceptable (x = 74 
for 51 dof) because the spectrum is more curved than the model, 
suggesting the presence of excess absorption. We then keep the 
Galactic column fixed, but add a second (neutral) absorption com- 
ponent. The statistical quality of the fit improves dramatically and 
wc obtain x 2 = 37 for 50 dof. We measure an excess absorption 
column density of 2.3± u ;| x 10 21 cm , while the power law has a 
slope T = 1.95t ;J|. The spectra, best-fit model and residuals are 
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shown in Fig. 3. We point out that the observed amount of excess 
absorption is typical for ULXs, suggesting that part of it originates 
in the ULX environment (Swartz et al. 2004; Winter, Mushotzky & 
Reynolds 2006). The X-ray spectral slope also compares very well 
with other ULXs and lies at the boundary between "low-state" and 
"high-state" sources (as classified by Winter et al. 2006). 

The observed 0.5-7 keV flux is 8.2 x 10~ 14 erg cm~ 2 s" 1 . 
By extrapolating the model up to 10 keV, the source has a 0.5- 
10 keV flux of 9.8 x 10~ 14 erg cm" 2 s" 1 . Note that the observed 
variability (see Fig. 2) implies that during the observation ULX 1 
reached a flux of about ~ 1.2 x 10~ 13 erg cm -2 s _1 in the 0.5- 
10 keV band. 



3.1 The extreme luminosity of ULX 1 

The excess absorption we detect and the lack of a bright opti- 
cal counterpart (see discussion below) makes it unlikely that the 
source is a foreground object, though we cannot exclude it. If so, 
the X-ray luminosity of the source can be computed as a lower 
limit under the hypothesis that it is at least at the distance of 
MCG-03-34-63 (z=0.0213). By using our best-fitting absorbed 
power-law model, it turns out that ULX 1 has an intrinsic lumi- 
nosity Lo.5-7 = 1.1 x 10 41 erg s~ 4 in the 0.5-7 keV band and 
Lo.5-10 = 1.3 x 10 41 erg s _1 once extrapolated to 10 keV (with 
H = 70 ks s~\ A = 0.73, and q = 0). If the light curve 
variability is taken into account, the source reached a luminosity of 
~ 1.35 x 10 41 erg s" 1 (~ 1.6 x 10 41 erg s" 1 ) in the 0.5-7 keV 
(0.5-10 keV) during the XMM-Newton observation. 

We point out that of the ~230 ULXs detected and catalogued 
up to now (e.g. Liu & Mirabel 2005; Ptak et al. 2006), only a few 
have luminosities exceeding 3 x 10 40 erg s" 1 . Luminosities of 
the order of 10 41 erg s _1 or above are extremely unusual, as also 
shown by the luminosity cutoff in the cumulative luminosity func- 
tions (Swartz et al. 2004; Gilfanov, Grimm & Sunyaev 2004). This 
makes ULX 1 a rare, if not unique, source among the other known 
ULXs. 



3.2 Adding a thermal component 

To test the similarity of the X-ray properties of ULX 1 with stellar- 
mass accreting black holes, we also tried two different models: a 
multi-colour disc (MCD) spectrum (DISKBB) plus power law, and 
a Comptonization model (COMPTT) that describes Compton up- 
scattering in a corona of optical depth r and electron temperature 
fcT c i from a Wien distribution of soft seed photons from the disc 
(with temperature kT). It is worth noting here that the detection of 
a thermal component could give important (though not conclusive) 
insights on the compact object mass since, for standard accretion 
disk, Tdisc oc (objects hosting black-holes of ~ 10 Mq 

have typical temperatures of 1 keV). A cool thermal component 
from the accretion disc may thus provide some indication for a par- 
ticularly massive black hole (see e.g. Miller, Fabian & Miller 2004). 

However, when applied to the pn and MOS data, neither 
model provides a statistically significant improvement with respect 
to the power law one. We can only measure 90 per cent upper 
limits for the temperature of the putative thermal component of 
kT < 0.70 keV (MCD model), and kT < 0.35 keV (COMPTT). 
Thus, the data are consistent with a cool disc, but the thermal com- 
ponent is by no means required by the data. 



4 AN EXTREME ULX OR A BACKGROUND AGN? 

The extreme luminosity estimated above assumes that ULX 1 lies 
within its apparent host galaxy MCG-03-34-63. Luminosities ap- 
proaching ~ 10 41 erg s -1 can also be reached by extremely power- 
ful supernovae occurring in dense environments, but given that their 
emission either fades or remains constant over a ~1 yr timescale 
(e.g. Schlegel 1995), the observed X-ray short-timescale variabil- 
ity rules out this hypothesis. However, ULX 1 could be a back- 
ground X-ray source. In this case, the inferred luminosity is only 
a lower limit and a background AGN could well be invoked to ex- 
plain the X-ray power of the source. This hypothesis can only be 
tested observationally by obtaining an optical spectrum of the coun- 
terpart (if any), but we can nevertheless provide two arguments in 
favour of an association with MCG-03-34-63. 

4.1 The log TV-log S estimator 

The number of expected X-ray sources within a given area can be 
estimated by assuming that our field is characterised by the same 
log AT-log S as observed from deep XMM-Newton exposures of 
the Lockman Hole (Hasinger et al. 2001). Our ULX candidate has a 
0.5-2 keV flux of 2.8 xlO" 14 erg cm" 2 s" 1 and the Lockman Hole 
log TAT-log S implies the presence of about 30 sources per square 
degree with the same or higher flux (we point out that 30 sources 
is a conservative estimate). Since the D25 ellipse of MCG-03-34- 
63 covers an area of 2.4 x 10~ 4 square degrees, the number of 
expected sources within D25 is 0.007. In the 2-10 keV band, we 
measure a flux of 7 x 10" 14 erg cm -2 s" 1 which again corresponds 
to about 30 sources per square degree providing the same small 
number of expected sources. 

As mentioned, this is only the number of expected sources 
within a given area (the D25 ellipse of MCG-03-34-63) and not di- 
rectly the probability that ULX 1 is a background AGN. We provide 
the above estimate for comparison with previous works on other 
ULXs, but we must point out that an accurate estimate of the true 
probability that an X-ray source within a galaxy is a background 
AGN has yet to be produced, and is well beyond the purpose of this 
Letter. 

4.2 The X-ray to optical flux ratio 

Another possible way to constrain the nature of the source is by 
considering the X-ray to optical flux ratio (X/0= -Fx/Jbpt). 
AGN have typical X/O below 10 and detailed studies of the op- 
tically faint sources in the XMM-Newton HELLAS2XMM and 
Chandra deep fields surveys reveal that less than 5 per cent of the 
hard X-ray selected sources have X/O higher than 90 (Mignoli et 
al. 2004; Civano, Comastri & Brusa 2005). 

We have computed the X/O for the brightest, and most likely 
optical counterpart, HST source in the XMM-Newton error-circle 
(i.e. the source within the yellow box in the bottom right panel of 
Fig. 1). From the HST image, we estimate a magnitude rriF606w = 
24.5 giving rise to R~ 23.8 (see Holtzman et al. 1995 for the 
conversion between HST and Johnson magnitudes). By using the 
definition log X/O = log Fx + 0.4 R + 5.61, we then estimate 
X/0~ 94 which represents a lower limit on the X/O of ULX 1 
(where Fx is the 2-10 keV flux). As mentioned above, sources with 
X/0> 90 represent less than 5 per cent of the hard X-ray sources 
population in Chandra and XMM-Newton deep fields. Moreover, 
the X-ray and near infrared data of high-X/O sources in the Chan- 
dra deep fields and the HELLAS2XMM survey identify most (if 
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not all) sources as heavily absorbed AGN (Civano, Comastri & 
Brusa 2005; Mignoli et al. 2004). In the case of ULX 1, the steep 
photon index and moderate absorption are inconsistent with this 
scenario. The detection of an unobscured background AGN with 
X/0> 94 would be a very rare event on its own, and its detec- 
tion within the D25 ellipse of a spiral galaxy seems highly unlikely, 
supporting the idea that ULX 1 is not a background AGN but more 
likely an X-ray binary. It should also be pointed out that the ULXs 
that have been later classified as background AGN had standard 
rather than high X/O (Foschini et al. 2002; Masetti et al. 2003). 



5 CONCLUSIONS 

If ULX 1 is indeed associated with its apparent host galaxy MCG- 
03-34-63, its luminosity greatly exceeds that of a neutron star, and 
we have to consider an accreting black hole. In this case, by con- 
sidering as a template stellar-mass accreting black holes in our 
Galaxy, the bolometric luminosity is likely to be a few times higher 
than the 2-10 keV luminosity. By extending our spectral model, 
we estimate an intrinsic luminosity of ~ 3 x 10 41 erg s _1 in the 
0.2-150 keV band, which we take as a reasonable proxy for the 
bolometric luminosity of ULX 1 . In the absence of beaming, the 
Eddington limit can be used to infer the presence of a > 2300 Mq 
black hole, i.e. ULX 1 would be powered by an accreting IMBH. 

To reduce the mass of the compact object within the accept- 
able range of stellar-mass black holes, the intrinsic luminosity has 
to be boosted by a factor > 115, which could be due to beaming. 
The simplest candidate mechanism to produce X-ray beaming in an 
accreting source is to consider a thick accretion disc with a smaller 
optical depth over a limited range of angles around the rotation axis 
producing thereby a central funnel. Such a geometry has been in- 
vestigated by Madau (1998) who found that if maximal beaming is 
considered, the effective observed luminosity (L c g ) can exceed the 
Eddington limit by a factor of ~23 (see also Misra & Sriram 2003 
who derived a much smaller factor of ~5). Thus, even by consid- 
ering maximal beaming with L c ff = 23 x Z/Edd, ULX 1 does still 
require a black hole with mass > 100 Mq. The mass of the central 
black hole could be in principle reduced to 1O-2OM0 by consid- 
ering maximal beaming and a source which is super-Eddington 
by a large factor (5-10), which seems a rather finely-tuned choice 
of the parameters. However, Begelman, King & Pringle (2006) re- 
cently pointed out that such a scenario is possible though it would 
predict a strong soft thermalized component (not required by the 
X-ray spectrum of ULX 1). Beaming can also occur if the X-ray 
emission is associated with a relativistic jet whose axis is aligned 
with the line of sight, as proposed by e.g. Reynolds et al. (1997) 
for ULXs. Indeed, a source with the jet axis aligned within a few 
degrees with our line of sight and with a sufficiently high Lorentz 
factor (7 > 4.3) could account for the observed luminosity with a 
sub-Eddington jet from a stellar-mass black hole. 
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